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ABSTRACT 

We discuss aspects of the low energy phenomenology of the MSSM, in the large 
tan P regime. We explore the regions of the parameter space where the ht and 
hb Yukawa couplings exhibit a fixed point structure, using previous analytic so¬ 
lutions for these couplings. Expressions for the parameters and Ab and the 
renormalised soft mass terms are also derived, making it possible to estimate ana¬ 
lytically the sparticle loop - corrections to the bottom mass, which are important 
in this limit. 
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1 Introduction 


Among the extensions of the Standard Model (SM), supersymmetry seems to provide 
the best grounds towards a unification of the fundamental interactions. This is not 
only because of the natural solution to the hierarchy problem |||], but also because 
of the correct prediction of sin^ 9\y , and the convergence of the three gauge couplings 
at a point, at a scale 0(10^®) GeV, unlike in the non-supersymmetric grand unified 
schemes Q. The simplest supersymmetric extension of the theory is the Minimal 
Supersymmetric Standard Model (MSSM) which has the minimal number of fields and 
Yukawa couplings that is consistent with the content of the SM. However, even in the 
MSSM several arbitrary parameters also exist. These are the initial conditions for the 
Yukawa couplings and the scalar masses, which are expected to be fixed at an even 
more fundamental level, like string theory. A simplified approach to reduce the number 
of these variables, would be to assume universality of all scalar masses at the unification 
scale. Then, one is left with only five new arbitrary parameters in addition to those 
of the non-supersymmetric standard model. However, while this is consistent with the 
low tan/3 regime, in the large tan/3 limit of the theory (where ht ~ hh), in order to 
get the correct radiative electroweak symmetry breaking pattern one is forced to 
depart from universality. This latter case, i.e. large tan/3 and non-universal boundary 
conditions for the scalars, appears naturally in many string derived models. In fact it is 
generally expected that irrespectively of the fundamental theory, the low energy model 
will look much like the supersymmetric standard model with non-universal boundary 
conditions for the soft terms. Further theoretical expectations suggest that the Yukawa 
couplings of the third generation are of the order of the common gauge coupling at the 
unification scale. Moreover, specific grand unified groups suggest equalities of the form 
ht = hb = hr at M[/. 

It is interesting to investigate whether radiative corrections can determine the Yukawa 
couplings and possibly the soft masses of the effective low energy theory. For example, 
the Yukawa couplings are running quantities from the unification scale Mjj down to low 
energy. If they are relatively large at Mjj, their low energy values exhibit a fixed point 
structure!^] and they are rather insensitive to their initial conditions. The analysis 
of the above becomes more interesting by the fact that there is recent experimental 
evidence for a top quark with a mass of O(180GeH) which is compatible with the 
prediction of the fixed point structure. 

In the present letter we investigate the fixed point structure of previous analytic solu¬ 
tions for the ht , hb couplings i§0- We further extend the existing analytic solutions 
1^ for the scalar masses in the case of tan/3 3> 1, by including the contribution from the 
At and Ab terms. This is of particular importance not only for the determination of the 
scalar mass spectrum itself, but for the correct theoretical computation of the fermion 
masses. In particular, when tan/3 S> 1 the bottom mass receives large corrections from 
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sparticle loops |H 10|. The tau lepton receives also corrections of the same type but 


they are less signihcant. A precise determination of these corrections would require 
the knowledge of the scalar masses involved in the computation. An analytic approach 
in both Yukawa and soft mass terms would make possible a systematic exploration of 
these corrections. 


2 Yukawa coupling fixed point 


In this section we are going to use the results of |8| for the top and bottom Yukawa 
couplings, in order to identify the regions where these couplings exhibit a fixed point 
behaviour and give simplified expressions that describe these fixed points. The renor¬ 
malisation group equations for the top-bottom system (when ignoring the h-j- Yukawa 
coupling), read 


hf = 


dt 

d ,2 


+hl - Gq'^hf 

^{hj + 6hl-Gn}hl 


where 


( 1 ) 

( 2 ) 

(3) 


Gq — Y dgSi , Gd — Y 

i=l i=l 

Here t = InQ, where Q is the energy scale, Cq = |^)3, and = |^)3, for 
C/(l), SU{2) and SU{3) respectively. Ignoring small U{1) corrections, this system can 
be solved. Dehning the parameters x, y through = 7 qX, hi = 'JqV, where 




®p{8?r°«<‘'K}=ii(^ 


one can make the transformation 

ko 


kt) 


u = 


j=i 


dl= ® 


Igdt 


(4) 


(5) 


(x — ’ Stt^ 

where u = x — y and I{t) = ^Q{t')dt'. The parameter 

depends on the initial conditions xq = hi q and yo = one forms a differen¬ 

tial equation for the new variable u, which can be solved in terms of hyper geometric 
functions. In particular, the t — b Yukawa coupling solutions can be expressed in terms 
of the variable u as follows Q : 


hi = 7t7qw(V 1 + -k 1) 


hi — -7i)u;(Vl + U — 1) 


(7) 

( 8 ) 


2 








Note that we have “restored” the symmetry between the differential equations 
and the solutions by replacing the gauge factor Jq with 7 I, in (^). 

It is interesting to search for particular combinations of the Yukawa couplings which 
are rather independent of their initial values. To start our investigation we first obtain 
a simplified formula for the function Lv(t) = x — y, in the case ^ 1 


uj{t) 


a^o - yo 


(9) 


where 2 F 1 is the value of the hyper geometric function at u = uq. Then, we can use the 
relation between x,y variables|^] 


1 

t- 

1 

( xy \ 

- yo^ 

^xoy^) 


to obtain the following expression, 


, -_ _ y/XQyo 

“ 2F^ + l^/xWoI{t) 


( 10 ) 


( 11 ) 


The above equation can also easily be obtained by a direct multiplication of the solutions 
(ID substituting u from (!)• 

Note first that for a small difference between ht and hh, 2 F 1 ~ 1, while the integral 
I > 10. Thus for ^x^yo = htfihbfi > 1 we can write 


hth 




( 12 ) 


This last expression tells us that we can get an approximate, model independent pre¬ 
diction for the product of ht,hh couplings at the low energy scale provided we start 
with relatively large and comparable htfl,hb^Q values at Mu. We note in passing, that 
it is possible to use Yukawa coupling constraints obtained at the Unification scale to 
eliminate one of the two parameters. Suitable constraints combined with ( 12 ), can 
determine the absolute values of both couplings. In |11] for example, it is shown that 
within spontaneously broken N = 1 supergravity models two generic types of such 
constraints are obtained. For a superpotential W = while assuming 

different scale structures for the various fields in the Kahler potential, there are mul¬ 
tiplicative (duality invariant) constraints of the form YliAisic diAiaic ~ course 

hiAigic 7 ^ 0. Even if this constraint were applicable for the two Yukawa couplings 
which interest us here, it could not be useful however, due to the fixed point property 


of the product in relation ( 12 ). A more interesting situation arises in cases where the 
constraint applies to the ratio of the Yukawa couplings. A simple example is shown in 
[ 111 for two couplings with a Kahler potential having the two fields in the same no-scale 
structure. Thus if the higher theory can give a prediction about the ratio of the two 
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Yukawas at Mjj, we can use this in conjunction with the analytic solution to extract 
information about the Yukawas at the weak scale. 

Then the low energy values can be obtained using the relation (0) and the low energy 
ratio obtained from the analytic solution in the limit htfihbfi > 1, 

§ » 2|^^±i (13) 

For the limit of interest, u is given by the approximate formula 

+ (14) 

with 

no = (2ro/(rg - 1)) (15) 

and ro,po the ratio and product of couplings respectively. Thus in the case 

of ro ^ 1, u ^ oo, and the ratio of the Yukawas runs also to a fixed point value 
determined approximately by the ratio ’JqIjd- 

In Figs, la, lb and Ic, we show the fixed point structure of the Yukawa couplings and 
their product, when we use the analytic expressions without making any approximation 
on the hyper-geometric function|^. We have taken a common gauge coupling ac = 
1/25.0 at a unification scale me = 1.35 • 10^® GeV and a supersymmetry breaking scale 
« 200 GeV, leading to an as(M^) 0.112. The couplings are presented at the top 
mass. 

Fig. la shows ht-low energy coupling versus the GUT /i/ values for the ratios h^/h^ = 
1.2 and 1.001 denoted with stars and crosses respectively. For a wide /i/ range, the 
estimated ht{mt) values differ at most by 1%. 

Fig. lb shows the bottom coupling versus h^. Here we took the same region of initial 
values for /i/ for each of the two cases quoted, thus the two lines are interrupted before 
touching the contour. In Fig. Ic we plot the product hbht for the same input ratios as 
in fig la. In Fig. Ic, a rather interesting fixed point property is exhibited if > 4 
were the estimated low energy values differ in less than 1%. On the other hand, it is 
remarkable that different initial /i/ /ratios accumulate exactly on the same curve at 
the weak scale. This is in a very good agreement with the equations that we derived 
for the description of the fixed points in terms of functions of the gauge couplings only. 

^The fact that we use expressions up to only one-loop, as well as the approximations made 
in order to derive the analytic formulas for the couplings, cause small errors in the numerical 
values of ht and hb which however do not alter the validity of the results. 
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3 Analytic Solutions for At and At 


The differential equations which govern the evolution of the soft scalar masses of the 
third generation and the two Higgs mass parameters niHi 2 ; known. In order 

to solve them, we need first a solution for the At^b,T trilinear mass terms. Ignoring the 
r-Yukawa coupling, we can write the At and A^ evolution equations as follows 


dAt _ 1 

dt Svr^ 

dAb _ 1 

dt Svr^ 


(^6htAt + h^Ab + GQmi/2j 

^6hfAb + hfAt + Gsmi/ 2 ) 


where Gq and Go are given by 




dvr 

ac 




2 = 1 


dvr 




“Gi.i 




(16) 


(17) 


where we have taken into account that the gluino masses are given by 

Mi Ri mii2— (18) 

OG 

og being the common coupling at unification scale. Since the small corrections that 
arise from the U{1) factors may be ignored we have Gq k, Gd. 

To solve this system we follow the lines of [^, where the system A4f/ — A4^ (ignoring At, 
Ab) which has a similar structure has been solved |^. We initially separate the running 
of the Yukawa couplings by rewriting At and Ab as 


At — rXt, Ab — aXb, 


where 


T = 




3 H 
to 


Defining the 2x2 matrix 

Tdit) = 7Q(o'r)T^ 

and the function 


0 2 / 0 ( 7)12 

^o(^)^ 0 


h{u) = 


xo \ 5 / _ 1 

2/0 / V V'i- + u + 1 


the At — Ab system can be written as 


A( 

du I Xb 


10 \/v? + u 


h{u) 


Ku, 




Xt \ ^ 

Xb I Svr^ du I Hb 


We come back to the system Aiff — Ml in section 4. 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 

( 23 ) 
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Note that h{u) in the large tan(3 case , (i.e. htfi ~ /ib,o) is approximately constant in 
most of the range of integration. 


Then, in the case n S> 1 we can write an approximate analytic solution as follows 


At 

Ab 


2 |(^? + hoA^)p + {A^ — — rni/2 < Jt > 

|(^t + ^oAb)p — {A^ — hoA^)-'^ — mi/2 < Ja > 


(24) 

(25) 


with 


, tand) . 1 
P=i - 


' tan(j)Q 


)6 sin2(j) = (1 + tt) 2 


oo 


while ho = h{u 

Furthermore, in the limit tt 3> 1, the integrals < Jr^a > are given by 

^‘°Gq(s) 1 


< Jt >= T 


-ds 


. (26) 

Jt t{s) p{s) 

and similarly for Jo-- A simple inspection of the above formulae shows that for rea¬ 
sonable initial At^b values the terms proportional to mi /2 dominate. Thus, to a good 
approximation we may write At ~ —?t7-i/2 < Jt > and Ab ~ —^- 1/2 < Ja >■ 


The semi-analytic expressions (24-^) are going to be used in the following sections, 
in order to compute the contributions to sparticle masses, as well as the corrections to 
the bottom mass from superparticle contributions. 


4 Predictions for sparticle masses and compar¬ 
ison with the exact solutions of the RGE 

Having obtained the solutions for At and Ab (thus for Aifj and we may calculate 
contributions to superparticle masses, by solving the system 

Mij = tuqt. '^'u + 

M\, = -|- rh\, + + A\ (27) 

with initial conditions AJfijD'jQ = i{u,D)'^‘o : where mo is a common scalar mass at the 
unification scale, = Chi + Taking into account the 

renormalisation group equations for the squark and Higgs fields^ , we obtain 

^ (28) 

^ In the renormalisation group equations of the squark and Higgs fields, for non-universal 
initial conditions, contributions proportional to where S{t) = ^ are obtained. 

However, these contributions cancel in the equations that describe the sums M'fj and Ai'jj, due 
to the invariance of the U{1) Yukawa Lagrangian. The non-universal initial conditions are still 
manifest, through the factors . 
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dM% 

dt 




( 29 ) 


where Gfj = Gq + Gh 2 + Gijc and G% = Gq + Ghi + Gb<^- Taking into acconnt the 
At and Af, contributions, the solution of the system to first order reads 


(30) 


Ml ^ 


ml 

2 1 

Ml ^ 

a 

ml 

2ho 


where +^Q + ^b<^- < U > is given by 


< It >= T 


Jt T{t') p{t') 


Finally, 

< Ar >= 2 t 

Here, At is approximately given by 


*0 / At(s) A[{s) 


t V t{s) p{s) 


ds 


and 


A,{t) ^ 




(32) 

ints r - 

-> a and Gu - 

-^Gd. 



(33) 

ds 


(34) 

Ggit) 

T{t) 

1 

pit) 

(35) 


where primes denote derivatives with respect to the logarithm of the scale t = log{p). 
Similarly, < > is found by the proper substitutions. 


Integrating the sums, we obtain 
y 2 a A'i 


M-u — -Mbq — Gu{t)m^i2 — —QJu — Jd + At — AtQ 
M\) — — CD{t)m\i2 = —Jjj ~ QJd + A^ — A^q 


(36) 

(37) 


with Ji{t) = / h‘]M.^jdt, I = U,D. Now, the unknown integrals Ji{t) can be expressed 
in terms of the already calculated functions, their initial conditions and known gauge 
functions, from the simple algebraic system (^,37). As an example, for the up-squark 
running mass squared we have 


— (?Q + C'Q(t)^i/ 2 )m,Q — Jnit) — Ju{t) + I's 


(38) 


with I'g representing the integral of the S'-contribution in the case of non - universality 
and 


CQ{t) = E 


3 Q 


^ 2hia. 


*G 


(39) 













"^gluino 

252 

252 

252 

252 

rriQL 

182 

185 

189 

188 

mtj, 

185 

188 

193 

191 


153 

157 

162 

160 

At 

143 

147 

- 

- 

Ab 

144 

146 

- 

- 


Table 1 : Some sparticle masses (in GeV) obtained using analytic (second column) 
or numerical (third column) solutions. The last two columns refer to the corresponding 
cases but when the terms are ignored. 

In table 1 we compare the superparticle masses obtained from the analytic formulae, 
with those when solving numerically the renormalisation group equations: In the first 
column, we give mass terms, for a numerical solution of the renormalisation group 
equations, and at a scale ~ 250 GeV, when the heavier superparticle, which in our case 
is the gluino, decouples from the spectrum. The initial conditions we take are 

A^t=Al = 1 , 

hO = 2.0, h^Jhl ^1.1, 

= 4.0, = 1-0) 

Crfc = 1.0, ^n- = 1.0, ^9 = 1.0, 

and ignore for simplicity the 1$ contribution to scalar masses (of course we do the 
same in the analytical solutions). In the second column appear the analytic solutions, 
when we include the At, Ah contributions, while in the third the solutions when these 
contributions are neglected. Finally, for comparison, in the fourth column we give the 
superparticle masses that are found numerically by the renormalisation group equa¬ 
tions, when we ignore the At^ contributions. Here we have taken niQ = mi /2 = 100 
GeV and the superparticle masses are given in GeV. The effect of the r coupling has 
been ignored in the numerical solutions of the renormalisation group equations as well, 
for a better comparison. Gomparing the analytic with the numerical solutions, we find 
that the total relative error is at most 2 — 3%. Part of this small error in the case where 
we include the At^b contributions arises because we wanted to keep the expressions as 
simple as possible and therefore have ignored the contribution of the homogeneous part 
of the At^b solutions in the superparticle masses (but not in the calculation of the H’s 
themselves). Had we included the effect of the r Yukawa coupling in the numerical 
solution of the renormalisation group equations, the shift to the results is at most 10%. 
These observations are in agreement with the ones obtained in the Appendix A of 0. 

Let us hnally make the following observation: In the equation for the product of the 
Yukawa couplings, we see that the larger the Yukawas and the smaller the difference 
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between ht — hb, the stronger the fixed point behaviour becomes. Moreover, in the 
equations for the scalar masses, we observe that the smaller the difference between 
— h^, the less important the contribution of the homogeneous part of the solution 
(which includes the main dependence from initial conditions for the Yukawa couplings, 
as the integrals < > undergo a smaller change) becomes. This indicates that for low 

energy Yukawas which correspond to and values close to the non-perturbative 

region for the couplings, one in principle expects to obtain the minimal variations of 
the effective potential under small displacements of the couplings. In 0 it has been 
shown that the effective potential exhibits a minimum related to the infrared fixed line, 
while for non-universal boundary conditions for the scalars, minima which correspond 
to ht = hb will also be favored. This we think can be understood by looking at the 
regions where the expressions we have derived exhibit the strongest attraction to the 
fixed points. 


5 Corrections to the bottom mass 


It has been found that in the large tan/3 regime, there are large contributions to the 
running bottom quark mass ruf,. These are given by 


m 


mb = hbVi (1 + 6mb). 


(40) 


where 6mb = ebtanf3 and 

and the integral function I(a,b,c) is given by 

a61n(a/6) + 6cln(6/c) + ocln(c/a) 


I (a, b, c) = 


(a — b){h — c){a — c) 


(41) 


(42) 


with Mg and {'niii) being the gluino and sbottom (stop) eigenstate masses respec¬ 
tively. 


To calculate these corrections we substitute the relevant expressions in the above formu¬ 
lae. In Q, expressions for the scalar masses have been derived (and these expressions 
are not sensitive to At, Ab). Here, we have derived an equation for At, while the gluino 
mass is approximately given by (|^). For the low energy value of /x we obtain the 
analytic expression 

~ (43) 

while the renormalisation group equation can be integrated to give 



( 44 ) 
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taniS 

mo = mi/2 

fi 

At 

/i(io-®) 

/2(10-®) 

Snib 

58.1 


124 

-147 

12.8 

18.1 

0.40 

59.3 


192 

-208 

6.2 

8.5 

0.41 

60.0 

■■ 

253 

-264 

3.8 

5.1 

0.42 

60.3 

250 

310 

-316 

2.6 

3.4 

0.43 

60.6 

300 

365 

-364 

1.9 

2.5 

0.43 


Table 2 : Bottom mass corrections for ho ~ 2.0. 


tan[3 

mo = mi/2 

/i 

At 

/i(io-6) 

/2(10-®) 


55.5 

HjjjH 

119 

-153 

12.1 

18.0 

0.34 

56.6 

msm 

184 

-216 

5.9 

8.4 

0.35 

57.1 


243 

-274 

3.5 

5.0 

0.36 

57.4 

250 

297 

-328 

2.4 

3.4 

0.37 

57.5 

300 

348 

-377 

1.8 

2.5 

0.37 


Table 3 : Bottom mass corrections for ho ~ 1.0. 


with {c^}i=i, 2,3 = {f,3,o|. 

Then, we find the magnitude of the bottom corrections, for solutions at the fixed point 
{ha = 2.0), as well as for small deviations from the fixed point {he = 1-0), while we 
keep h^/h^ ~ 1.1. The relevant quantities appear in tables 2 and 3. The masses in 
the tables are given in GeV. 

As we increase the supersymmetry breaking scale, tan/3 slightly increases, in order 
to get the same low energy parameters. (At the same time the uniheation scale 
drops slightly, while the inverse gauge coupling at the uniheation scale increases, 
by a small amount). 


6 Conclusions 

In this letter we have used simplihed analytic solutions for the ht, hb Yukawa cou¬ 
plings in order to study the MSSM in the large tanl3 regime. We have explored the 
regions of the parameter space which lead to a hxed point structure and derived 
the evolution of the Yukawas towards these hxed points. Using this information, 
one may identify the regions for the initial values of the Yukawa couplings which 
lead to the strongest attraction towards these infrared hxed points. Under these 
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considerations, top-bottom Yukawa coupling equality, and values of the couplings 
close to the non-perturbative regime seem to be favoured. Finally, we obtained 
corrections on the renormalised soft mass terms due to the evolution of the tri- 
linear parameters At and Ai,. Using these results, we estimated analytically the 
sparticle loop - corrections to the bottom mass, which are important in the large 
tan/5 scenario. In agreement with previous calculations we hnd that the maximal 
corrections arise at the hxed point. 


G.K.L. would like to thank the group of Centre de Physique Theorique de 1’ Ecole 
Polytechnique for kind hospitality. 
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